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Abstract 

It is demonstrated that the zero frequency conductivity, or Drude weight, 
can be expressed in terms of a geometric phase. The resulting expression is 
then cast in terms of ground state expectation values over the total momen- 
tum shift and total position operators, leading to easy numerical access. The 
criterion for insulation is demonstrated to be the localization in the many- 
body configuration space relative to the localization of individual charge 
carriers. 

The physical significance of geometric phases resulting from adiabatic cycles 
was first emphasized by Pancharatnam [[T| and Berry [2]. Since the appearance of 
these seminal papers the Berry phase has become ubiquitous in describing essen- 
tial physical phenomena (31 IH ranging in applications from quantum chemistry 
to topological insulators. One important quantity described by a Berry phase is 
polarization in crystalline systems [|5]|6). In the modern theory polarization is cast 
in terms of the integrated current between two states connected by an adiabatic 
path. In numerical applications the polarization is easiest to calculate in terms of 
the ground state expectation value of the total momentum shift operator [0 [S]. 
This development has simplified the calculation of the polarization considerably, 
ans is now in widespread use in density functional calculations. 

An analogous formalism with appropriate modifications is also valid for the 
linear current response 0. Polarization and current are described in terms of a 
particular type of geometric phase, sometimes referred to as a Zak phase [10]. For 
these phases the adiabatic cycle is a result of the periodicity of the relevant pa- 
rameter space (the crystal momentum in the case of polarization, and the average 
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position in the case of the current). The phase in this case is a line integral over 
an open curve, from one boundary of the periodic parameter space to the other. 

The Berry phase in general cannot be cast in terms of a line integral over an 
open curve. The integrals usually enclose surfaces in two or more dimensions. 
A well-known example is the Hall conductance IfTTTl . where the relevant phase 
is an integral around the two-dimensional Brillouin zone. The phase associated 
with Hall conductance was the first physical example of a topological invariant 
(TKNN invariant IfTTTl '). A number of important examples exist of Hamiltonians 
tuned to produce finite TKNN invariants. The Haldane model lfT2l is a modified 
tight-binding model on a hexagonal lattice which exhibits quantized Hall conduc- 
tance, and was the first example which does this in the absence of an external mag- 
netic field. Other extensions of the tight-binding hexagonal model, the Kane-Mele 
model IPT51 FBll . exhibit different topological invariants, and are the fundamental 
models for topological insulation in two-dimensions. Up to now these types of 
phases have not been expressed in terms of a ground state expectation value as 
was done by Resta (71[8|| for polarization. 

In this paper we first derive a general procedure for evaluating the spread 
in quantities which correspond to phases of the Zak type. In the course of this 
demonstration we argue that the quantities which define the Berry phase can be 
interpreted as operators in the space of external parameters which describe a sys- 
tem. For example the Hall conductance corresponds to a commutator between 
such operators. Subsequently we show that the zero frequency conductivity can 
also be cast in terms of a geometric phase, similar to the one corresponding to the 
Hall conductivity, as a commutator of heuristic operators. The expression for the 
phase contribution is then cast in terms of the ground state expectation value of 
shift operators, facilitating numerical access. The connection between localization 
and insulation, first suggested in a seminal paper by Kohn lfT51l . is made precise: 
the Drude weight, the criterion for insulation, depends on the localization in the 
many-body configuration space relative to the localization of individual charge 
carriers. 

The most general way to obtain the Berry phase is to write it in the discrete 
representation, and then take the continuous limit. One considers a Hamiltonian 
which depends smoothly on some parameter In this case the ground state wave- 
function will also be a smooth function of the parameters For a set of points {£,•} 
one can form the quantity 

N-l 

= -Imln] _ [<T(^)|^ /+1 )), (1) 

7=0 

where ^(^n)) = |*F(£o)) (cyclic process). The phase (p is a quantity physically 
well-defined, since the arbitrary phases of the different wavefunctions cancel. We 
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can consider a closed curve in the parameter space £ and take the set of points 
{£/} to be along the curve in the order of the indices /. We can then take N — » oo 
corresponding to the continuous limit. Expansion in the distances between points 
results in 

<f> = i£ d£. OP(0|V f PF(#. (2) 



The line integral in Eq. © is similar to an expectation value over the oper- 
ator Nf. In fact in the modern theory of polarization [13, in which polarization 
is expressed in terms of a Berry phase, the quantity idg (with K the crystal mo- 
mentum) plays the role of a "heuristic total position operator" lfT6l . Similarly, the 
current can be expressed as a Berry phase [J9] consisting of a line integral over the 
operator id x (where X denotes the sum over positions of particles). 

If we were to consider the real part of the product of scalar products in Eq. 
dH) the continuous limit would be ill-defined. Indeed the first order term in the 
expansion (corresponding to the Berry phase) would disappear, and in the second 
order term the quantity d£ would be squared. It turns out that one can nevertheless 
obtain a physically well-defined quantity from expanding Eq. (OQ) up to second 
order. To show this we first consider the Berry phase corresponding to the polar- 
ization as written by Resta B71|8|. We stress though that the derivation is general. 
Resta showed that the average position over some wavefunction I 1 ?) of a system 
with unit cell dimension L can be written as 

(X) = --±-Im\n{y\e- iAK *\y), (3) 
AK 

where AK = 2n/(NkL), denotes an integer, X = is me sum °f me posi- 
tions of all particles. The spread (cr| = (X 2 ) - (X) 2 in position can be written 

o\ = -^Reln<¥k-'' Arf m, (4) 

The operator e' AKX is the total momentum shift operator which, as has been shown 
elsewhere [[T7Tl has the property that for a state ^(K)) with particular crystal mo- 
mentum K 

e- iAKjt MK)) = \W(K + AK)), (5) 

in other words it shifts the crystal momentum by AK. To use the shift operator we 
first write 

a% = ^Re ln<¥|e-' A ^|¥) Wt . (6) 

x N k AK 2 

We associate the state I^P) with a particular crystal momentum K , 

m = ito>)>. (7) 
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Using the total momentum shift the scalar product can be rewritten as 



(y(K )\e- iAKX mK )) = OFTOIWi)) 

= mKj)MK I+1 )), (8) 

where K I+i = K[ + AK. To show the last equation one applies the Hermitian 
conjugate of the total momentum shift to (^(Kq)] I times and the total momentum 
shift operator to \*¥(Ko)) I + I times and forms the scalar product. Thus we can 
also write 

N k -\ 

(y(K )\e- iAKjt mK )) N * = YlrnKtWiK^)). (9) 

/=o 

The points Kj form an evenly spaced grid with spacing AK in the Brillouin zone. 
Using this result the spread can be rewritten as 

2 Nk 

o\ = -^gRelnm^)|W +1 )>, (10) 
We now expand the scalar products up to second order as 

(v(KjW(k m )) = i + AKmKjmmKj)) + ^(^kwmkj)). (ii) 

Subsequent expansion of the logarithm up to first order leads to a first and second 
order term in AK. The first order term is given by 

N L 2 Nk ~ i 

-^Re V AKWmidMKj)). (12) 
2n to 

In the continuum limit (Nk —* 00 ) the sum turns into the integral which gives 
the standard Berry phase, but since this integral is purely imaginary it will not 
contribute to the spread. The final result for the spread is 

a\ = Yj A£<W/)l4PW)> (13) 



l r 

2?r J_ ff 



1=0 

•7T/L 

,2 



dK(¥(E)\dimK)). 

-n/L 

Eq. (TT3l) supports the idea that the quantity id K behaves like an operator. Indeed 
the line integral over the square of id K results in the second moment of the total 
position (the square of the polarization centered around the average polarization). 
We define 

1 = id K . (14) 
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The derivation leading to Eq. (TT3T > holds for phases of the Zak type, in which the 
cyclic nature of the adiabatic path results from the periodicity of the parameter 
space. 

Another example for which the above derivation is valid is the current resulting 
from a Peierls type perturbation. We consider a system with Hamiltonian //(O), 
with <D denoting the Peierls phase corresponding to the perturbation. £(0) = 
( V P(0)|//(0)| V I , (<1))) denotes the ground state energy. To keep the discussion as 
general as possible we first expand the Hamiltonian up to second order in <J> and 
obtain ^ 

H(0) = H(0) + <MT(0) + y//"(0). (15) 
For both continuous and lattice Hamiltonians the following identities hold: 

H'(0) = i[H,X] = J, (16) 

and 

H"(0) = i[J,X]. (17) 

where X denotes the sum over all positions, and J denotes the current operator for 
the sytem considered, both independent of the boundary conditions. The operators 
X and J are, of course, different for continuous or lattice systems. 

The total current is given by the derivative of the ground state energy with 
respect to the Peierls phase. Using the second order expansion of the Hamiltonian 
one can show that 

J ^ = = a ° + (¥(W(0)rF(4»>, (18) 

where 

a = ;<Y(Q)|[/,Z]PF(0)>. (19) 

It has been shown elsewhere [9| that the first term can be converted into a 
geometric phase. The result of that derivation is 

/(0) = aO-i- f dX(¥(®,X)\d x ¥(<I>,X)), (20) 
L Jo 

where ¥($, X) denotes the ground state wavefunction under a perturbing field O 
and with total position X. Wavefunctions with any X can be obtained by applying 
the position shift operator to some starting wavefunction with some total position 
X. One can write the discrete analog of 7(0) as 

7(0) = aQ> - -3-Im ln(¥(0, X Q )\e~ iAxk M^, X )), (21) 
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where X is the total position of the ground state wavefunction ^(O, X ) and AX = 
L/N x , with N x integer. The operator e~ lAXK is the total position shift operator HTS1 
13, and has the property that 

e ihxt |Y(X)> = MX + AX)). (22) 

Without loss of generality we restrict to an unperturbed system O = 0. For such a 
system the spread in current (cr 2 = (J 2 ) - (J) 2 ) may be written as 

cr) = ^Re \n(y(X )\e- iAxk I TO,)). (23) 
The continuous limit (N x — > oo) in this case is 



(J 2 ) - (J) 2 = -- I dX<T(X)|^| W)>, (24) 



^ Jo 



where X denotes the average position within the periodic cell. The expression idx 
can be interpreted as a heuristic current operator. We define 

J = id x . (25) 



We now consider the DC conductivity 11151 and show that it can also be ex- 
pressed in terms of a geometric phase. The DC conductivity is given by the ex- 
pression Q3] 



~ 7T 

D c = - 
L 



d z E{Q>) 
<9$ 2 



(26) 

$=o 



To obtain the Drude weight we take the derivative of the current 7(<B) with respect 
to O and obtain 

<)-tm = i r ^ (<W(0fJ Q| axY(<l>> ^ > _ <axT(0> ^| W0>X)>) _ (27) 



SO 2 L Jo 

The derivative d<j> can be set equal to the derivative in the crystal momentum d K , 
and subsequently an average over the Brillouin zone can be performed, which 
results in 

d 2 £(0>) 

ser ma+1 > (28) 

with 

y = — f f dtfdX«d* W W X)) - <fl z Y(ff, 
2^ J-^/l Jo 

(29) 



6 



One can see from Eq. (1281) that the Drude weight consists of the expectation 
values of two commutators; one of the standard current and position operators, 
the other between the total current and total position operators (Eqs. (fl4l) and 
(l25l)). The latter result are the contribution due to the boundaries, whose role in 
conduction vs. insulation was first emphasized by Kohn IH31 . The second term 
in Eq. (1281) corresponds to a geometric phase and it can be converted into a line 
integral around a closed curve by applying the Stokes theorem. 

We can also convert the Drude weight into an expression of the type derived 
for the polarization by Resta [|71[8]|. After an integration by parts and using the fact 
that the representation of the operator idx in reciprocal space is K we can write 

y = -- [' dKReK[mK)\d K ¥(K))l (30) 
Using a discretization of A^f = (2n)/(N k L) we can write the integral as 

7 = -TT X IRe[AK(W(K)\d K y(K))]. (31) 

For small AK we can replace the summand by 

T = -tt J] IR^V + AK(y(K)\d K *¥(K))]. (32) 

Nk I=-N k /2 

which is equal to 

2 N k /2-l 

T = -— 2] RzMViKWiK + AK)) 1 . (33) 

Nk I=-N k /2 

Using Eq. © and evaluating the sum on / results in 

y = ^Reln( v F(^)|e" ! ' A ^|T'(^)). (34) 

In simplifying the Drude weight we can go one step further. The quantity a is 
an expectation value over the commutator i[J,X\. The operators J and X are 
sums over the current and position operators of individual particles, as a result the 

N A 

commutator itself becomes a sum over individual contributions i Z ; ' Xj], with 
N p being the number of particles. The position of an individual particle can be 
represented as id~ kj , hence a takes the form 

a = iY j (V\[k j ,d kj ]m, (35) 
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which can, via the same reasoning used to arrive at Eq. (134T > can be converted into 

1 N " 

a = --Re In n<¥|<r^pp> (36) 

7=1 

in other words a sum over single particle contributions similar in form to y. The 
Drude weight then becomes 

m>\e- iAK *m 

D c = Re In — . (37) 

This results refine the tenet of Kohn [15] regarding the connection between insu- 
lation and localization. A DC conductivity of zero will result if the localization of 
the system as a whole is equal to that of the product of individual charge carriers. 

The above formalism can be applied, with minor modifications, to the Hall 
conductance, which is given by IfTTTl 

ie 2 C 

aH= 2^J dKxdK y ({ d ^ d ^> ~ <<VIV» (38) 

Here the expectation value is over the commutator expression i[d Kx ,d K ]. For a 
vector potential of the form (0, eBx) the Hamiltonian will have the form 

2m 2m 

Shifting the momenta in the y direction amounts to shifting the positions in the x 
direction, apart from the modification of the potential V(x, y) whose x component 
is also modified, but since that amounts to shifting the whole system the Hall 
conductance will be unaffected. One can then substitute the momentum shift in 
the y direction with a total position shift on the x direction and proceed as above. 

In this paper the DC conductivity (Drude weight) was expressed in terms of a 
geometric phase similar of the same form as the TKNN invariant. The expression 
was then converted into an expectation value over the ground state wavefunction. 
The Drude weight was thereby demonstrated to depend on the ratio between the 
localization of individual particles and that of the system as a whole. 

The author acknowledges a grant from the Turkish agency for basic research 
(TUBITAK, grant no. 112T176). 
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